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We study the possibility to test the Type I seesaw mechanism for neutrino masses at the CERN 
Large Hadron Collider. The inclusion of three generations of right-handed neutrinos (Ni) provides 
an attractive option of gauging the B — L accidental symmetry in the Standard Model (as well 
as an extended symmetry X = Y — 5(B — L)/A). The production mechanisms for the right- 
handed neutrinos through the Z' gauge boson in the U(1)b-l an d U(l)x extensions of the Standard 
Model are studied. We discuss the flavor combinations of the charged leptons from the decays of 
Ni in the AL = 2 channels. We find that the clean channels with dilepton plus jets and possible 
secondary vertices of the N decay could provide conclusive signals at the LHC in connection with 
the hierarchical pattern of the light neutrino masses and mixing properties within the Type I seesaw 
mechanism. 



I. INTRODUCTION 



The small but non-zero neutrino masses lead to a deep conjecture: Majorana nature of the neutrino 
masses may hold the key for a fundamentally different mass generation mechanism, although Dirac masses 
can be generated via the Higgs mechanism by introducing right-handed neutrinos with arbitrarily small 
Yukawa couplings. There are three simple scenarios where one can generate Majorana masses of the neu- 
trinos with renormalizable operators at tree level, and in agreement with the observations, the Type I |l|], 
Type II I2I, and Type III iQ] seesaw mechanisms. See also Refs. |4] and |5] for the simplest neutrino mass 
generation mechanisms using radiative corrections. 

Perhaps the simplest and best-studied mechanism for neutrino masses is the Type I seesaw, where one 
introduces at least two right-handed neutrinos (N). Adding in the corresponding large Majorana mass terms 
(M), one results in at least two light Majorana neutrinos with masses given as m? D /M. It is important to 
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mention that the inclusion of three right-handed neutrinos also provides an anomaly-free formulation for a 
gauged U(l ) B -L ffl- 

The non-ambiguous test of the Majorana nature of the neutrinos, and thus a possible test of the seesaw 
mechanism, will be the observation of the lepton number violation processes. The neutrinoless double beta 
decay is also a crucial test and one of the most sensitive probes. Since the CERN Large Hadron Collider 
(LHC) is going to lead us to a new energy frontier, searching for the heavy Majorana neutrinos at the LHC 
appears to be very appealing [7, ^]. However, due to the rather small mixing between the heavy neutrinos 
and the Standard Model (SM) leptons in a minimal Type I scheme, typically of the order | V^v| 2 ~ m u /M^, 
the predicted effects of lepton number violation are unlikely to be observable. On the other hand, if there 
are other particles beyond the SM that can mediate new interactions between them, the effects may be 
significantly enhanced. For instance, with the new gauge interaction U(1)b-l> the gauge boson Zb-l can 
be produced copiously at the LHC via its gauge interactions with the quarks. Its subsequent decay to a pair 
of heavy Majorana neutrinos may lead to a large sample of events without involving the small mixing angle 
suppression of N E,H]. The AL = 2 signals will directly test its Majorana nature; and the lepton flavor 
combination could probe the properties of the light neutrino mass spectrum and mixing pattern. 

In this paper, we investigate the possibility to test the Type I seesaw mechanism at the LHC in the context 
of two simple extensions of the Standard Model where one has an extra Abelian gauge symmetry. We focus 
our attention on scenarios with a U{\)b-l or U{l)x (X = Y — 5(B — L)/A), where B, L and Y stand 
for Baryon number, Lepton number and weak hypercharge, respectively. In order to cancel the anomalies, 
one just need to introduce three right-handed neutrinos, which are the source for the Majorana masses. In 
both scenarios one has a new neutral gauge boson, Z', which couples to the right-handed neutrinos. Then, 
one can expect large number events for the lepton number violating events due to the production and decays 
of the TeV Majorana neutrinos. The predictions of the heavy neutrino decays in each neutrino spectrum, 
Normal Hierarchy (NH), Inverted Hierarchy (IH) or Quasi-Degenerate (QD), are investigated in great detail. 
We find encouraging results for the LHC signatures to learn about the light neutrino properties. 

This work is organized as follows: In Section II we discuss the constraints on the mass and mixing 
parameters in the Type I seesaw mechanism from the current neutrino oscillation data. The predictions for 
the decays of the heavy neutrinos in the different neutrino spectra are presented in Section III. In Section 
IV we discuss the possibility to test Type I seesaw at the LHC through the same-sign dilepton channels. 
We summarize our findings in Section V. The mixing between light and heavy neutrinos are discussed in 
Appendix A. We provide the explicit expressions for these mixings in Appendix B. The minimal extensions 
of the Standard Model to U(1)b-l an d U(l)x are discussed in Appendix C. 
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H. TYPE I SEESAW MECHANISM AND PARAMETER CONSTRAINTS 

In the case of the Type I seesaw mechanism for neutrino masses one introduces at least two SM singlets, 
right-handed neutrinos, vr ~ (1, 1, 0), in order to generate two non-zero neutrino masses. In this case the 
relevant Yukawa interaction and the Majorana mass term are given by 

- Ci = Y U D l L Hu R + ^-ulCv R + h.c. (1) 

Here H = io2H* and the lepton number is broken in two units due to the presence of both terms. Now, 
integrating out the right-handed neutrinos one finds that the mass matrix for the light neutrinos is given by 

M v = m D M^ 1 ml, (2) 

where m/j = YJ^vq/V^, is the Dirac mass term and vq is the Higgs vacuum expectation value. Therefore, 
in this framework one could understand the smallness of neutrino masses, since the mass scale in the 
above equation could be large, M/v 3> YJ^vq, This is the so-called canonical Type I seesaw mechanism 
The mass matrix for neutrinos is diagonalized by unitary rotations as detailed in Appendix A. The three 
light neutrino masses can be expressed in the following way 

m = V PMNS M " V PAINS> (3) 

where m = diag(m\ , mi, m^) and Vpmns can be taken as the leptonic mixing matrix for the three gener- 
ation of light neutrinos jllll without the loss of generality. 1 Working in the basis where the heavy neutrino 
mass matrix is diagonal and using the Casas-Ibarra parametrization IU2I1 one can write mjj satisfying Eq. (f2]) 
as 



m D 



VpMNsm^VM 1 ' 2 , (4) 



where M = diag{M\, M2, M3) for heavy neutrino masses, and is a complex matrix which satisfies the 
orthogonality condition Q, T Q, = 1. It is shown in Appendix A that using the seesaw formula and the relation 
between the leptonic mixing one can find a formal solution for the mixing between the SM charged leptons 
(I = e, n, t) and heavy neutrinos (N = 1,2,3): 

V eN = Vpmns m 1 ' 2 n M- 1 ' 2 . (5) 

Therefore, for a given form of Q,, one can establish the connection between the heavy neutrino decays and 
the properties of the light neutrinos. The impact of the existence of the SI matrix on the decays of heavy 



1 The 3x3 rotational matrix is not exactly unitary when there are extra Majorana neutrinos, but it is a good approximation to 
equal it to the traditional Vpmns, see the formalism in the appendix. 
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neutrinos has not been studied before in collider phenomenology. Unfortunately, since the explicit form 
of this matrix is unknown one cannot predict the decay pattern of the heavy neutrinos with respect to the 
spectrum for light neutrinos. We will present a few well-motivated typical cases where one can hope to see 
the connection in each spectrum for light neutrinos. It is important, however, to realize that an underlying 
theory would pick only one specific form of tt. This (yet unknown) form would have definite prediction for 
the N decay patterns, through which the underlying theory could be revealed. 



A. Constraints on the Physical Parameters 
1. Neutrino Masses and Mixings 

In order to understand the constraints coming from neutrino physics let us discuss the relation between 
the neutrino masses and mixing. The leptonic mixing matrix is given by 

/ 



VpMNS 



C12C13 



C13S12 



e- iS s l3 \ 



-Cl2Sl3«23e 



C23S12 C12C23 - e S12S13S23 C13S23 



diag(e i$1 / 2 ,l,e^ 2 / 2 ) (6) 



\ S12S23 - e 1<5 Ci 2 C23Sl3 -C23«12Sl3e 1<5 - C U S23 C13C23 ) 



where = sin 0y, Cij = cos 6^, < 0y < ir/2 and < 5 < 2ir. The phase 5 is the Dirac CP phase, and 



<J?, are the Ma 



at 2d level 111 311 . are 



orana phases. The experimental constraints on the neutrino masses and mixing parameters, 
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(7) 
(8) 
(9) 
(10) 
(11) 



and m i < 1-2 eV. For a complete discussion of these constraints see reference 11411 . Following the 
convention, we denote the case Am^ > as the normal hierarchy (NH); Am^ < the inverted hierarchy 
(IH), and the quasi-degenerate (QD) spectrum where the lightest neutrino mass is larger than 5 x 10 -2 
eV. Using the above experimental constraints, one can expect to explore the allowed values for the Vin 
couplings and the heavy masses. From Eq. ©, we can obtain the general expressions of ^Ar(^2v) 2 tnat 
are collected in Appendix B. 
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FIG. 1: J^n \V(n\ 2 Mn /100 GcV versus the lightest neutrino mass for NH (left) and IH (right) in Case I (degenerate 
TV), assuming vanishing Majorana phases. 

2. Case I: Degenerate Heavy Neutrinos 

We firstly study the simplest case where the three heavy neutrinos are degenerate. This is a highly 



16b at the 



motivated scenario since it is strongly favored to generate successful resonant lepto-genesis |[15 , 
low scale. Using Eq. (IA15I ) and assuming degenerate heavy neutrinos we obtain the relation 

M £ (V; N ) 2 = {V^ MNS mVl MNS )u = {M v ) u , (£ = e^,r). (12) 

^=1,2,3 

We see that one can obtain simple relations for the heavy neutrino mixings and masses in terms of the light 
neutrino mass matrix independent of the unknown matrix f2, which in turn is given by the parameters from 
the neutrino oscillation data. One can thus predict the decays of the heavy neutrinos in each spectrum for 
light neutrinos. Note that in this degenerate scenario, we are unable to convert the constraints of Eq. (fl2l to 
predict ^T, N | V^/v| 2 in general. We can predict the decays of heavy neutrinos in terms of the other oscillation 
parameters only when all phases vanish since in this case the modulo square of the mixings (which govern 
the decay rate) are equal to the square of mixings (the left-handed side of Eq. (fT2l ). 



In Ref. 11711 . we have shown that using the experimental constraints on the neutrino mass parameters the 
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elements of the neutrino mass matrix has the following properties: 

Ml e <C M^,M T V T for NH, 

Mt e >M^,Ml T for IH, (13) 
M ee _ M w _ M rr for q D (14) 

Following the same approach, we plot the allowed values for the normalized couplings of each lepton flavor 
in this scenario in Fig. [H as a function of the lightest neutrino mass in each spectrum, the normal hierarchy 
(left panel) and the inverted hierarchy (right panel), assuming vanishing Majorana phases. We see two 
distinctive regions in terms of the lightest neutrino mass as expected. In the case m 1 ( 3 ) < 5 x 10~ 2 eV, we 
see the characteristic features 

E I^| 2 «E l^vl'E \ V tn\ 2 for NH, 

N N N 

E lW > E I WE l^l 2 fOT IH - 
N N N 

On the other hand, for mu^ > 5 x 10 -2 eV, the light neutrino masses enter the QD spectrum that leads to 

EiW«EiW«EiW- 

N N N 

Under this mass degenerate assumption, the mixing between the heavy neutrinos and the SM charged lep- 
tons simply reflect the features of the light neutrino mass matrix in the flavor basis, as seen in Eq. (PT2l ). This 
is an important model-prediction. It is important to emphasize that the results shown in Fig.[T]may be used 
to learn about the neutrino spectrum. 

3. Case II: Non-degenerate Heavy Neutrinos 

If we relax the assumption that heavy Majorana neutrinos are nearly degenerate in mass, then the com- 
plication due to the unknown matrix Q arises. The explicit parameterization of SI is presented in Appendix 
A, and the general expression for the relations among the parameters are given in Appendix B. 

For the purpose of illustration, let us take Q to be a real matrix. We could gain a general sense 
for the mixing parameters by varying the matrix elements of Q in the range of —1 < Wij < 1. 
We show |Vh| 2 Mi/100 GeV in this case in Fig. The predictions of \Ve 2 \ 2 M 2 /W0 GeV and 
| V£s\ 2 Ms/l00 GeV are almost the same. As one can see that qualitative features for both cases of 
NH and IH closely resemble that in Fig. Q] This is quite encouraging since the random selection 
of the model parameters do not seem to totally wash out the predicted features. To further explore 
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FIG. 2: |V^i| 2 Mi/100 GcV versus the lightest neutrino mass for NH (left) and IH (right) with Q = R12R13R23 and 
random matrix elements — 1 < Wij < 1, assuming vanishing Majorana phases. 

the model implications, we must choose a specific form of the O matrix, which should correspond 
to a particular theoretical incarnation in the right-handed neutrino sector. However, a large Majorana 



phase could alter the predictions Ill7ll in general. We will check on this point in the next section. 

Case Ila: tt = I 

In this simple scenario, we easily obtain transparent relations for the Nx mixings: 

|V e i| 2 Mi = micf 2 c? 3 « mic? 2 , (15) 

|V^i| 2 Mi = mi|si 2 C23 + ci 2 si 3 s 23 e l5 \ 2 ~ misl 2 C2 3 , (16) 

|y r i| 2 Mi = mi|si 2 S23 - c 12 s 13 c 23 e tS \ 2 w misf 2 *23' (I 7 ) 

and therefore 1 1 2 > |V^i| 2 , |K-i| 2 . In the case of /V2 mixing: 

|V r e2 | 2 M 2 = m 2 c 2 3 s 2 2 « m 2 s 2 2 , (18) 

|^ M 2| 2 M 2 = m 2 \c 12 023 - si 2 si 3 S23e t5 \ 2 « m 2 c 2 2 C23, (I 9 ) 

|V A r2 | 2 M 2 = m 2 |si2Si3C23e J ' 5 + C12S23I 2 ~ "T-2C 2 2s| 3 i ( 2 °) 
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and \ V e 2 




V T 2\ 2 - As for the N3 mixing, 



|^ e3 | 2 M 3 
|^ 3 | 2 M 3 
|K 3 | 2 M 3 




2 2 
^3Ci 3 C23 



™ 3 C2 3 . 



(21) 



(22) 



(23) 



and one can see |V^ 3 | 2 , \V T %\ 2 > \V e s\ 2 . These features are shown in Fig. [3] 

A few remarks are in order. First of all, with this choice of a diagonal matrix Q, the mixing angle squared 
| Vu\ 2 for Ni is always proportional to the corresponding light neutrino mass m,. Consequently, the relative 
fractions of the mixing to different lepton flavors are universal for both NH and IH. Secondly, the Majorana 
phases do not appeal - in | V£n\ 2 due to the special structure of fl. Thirdly, as seen in Fig. [3l the relative 
strength of the mixing to different lepton flavors for each N closely follow that for the light neutrino mass 
eigenstates. In fact, very much like the light neutrino mass eigenstate labeling, this should be the defining 
feature to label N±, N2 and iV 3 , if we do not like the less illuminating ordering Mi < Mi < M 3 . 



We choose to study yet another simple, but different form of the matrix, namely, with $7 as an off- 
diagonal unity matrix. As can be shown explicitly and supported by Fig. [4] the mixing features of | Vn \ 2 
and I Vts | 2 switch places with each other in both NH and IH, while | V12 \ 2 remain the same as in Q, = I case. 
If we recall the convention for the N labeling, this case is indistinguishable from Case Ha. In this case 
I V^jv | 2 are a l so independent of Majorana phases. This similarity can be generalized to a matrix of f2 which 
has only unity as entries. We would expect that the real situation could be a well-defined superposition of 
the three vertical panels as long as Q is real. 



The leading decay channels for the heavy neutrinos include Ni — ► ejW T , N — > v^Z and Ni — > 
Vjh(H). The amplitude for the two first channels are proportional to the mixing between the leptons and 
heavy neutrinos given in Eq. ©, while the last one is proportional to the Dirac-like Yukawa terms given in 
Eq. ©. 



Case lib: 17 = I g 



III. HEAVY NEUTRINO DECAYS AND LIGHT NEUTRINO SPECTRA 
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FIG. 3: | Vu\ 2 Mi/ 100 GeV, i = 1,2,3 versus the lightest neutrino mass for NH (left panels) and IH (right panels) in 
Case Ha (Q = I). 
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FIG. 4: \Vii\ 2 Mi/100 GeV, i = 1,2,3 versus the lightest neutrino mass for NH (left panels) and IH (right panels) in 
Case lib (0 = I oS ). 
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A. Decay Modes of Heavy Majorana Neutrinos with mass: Mi > M w 

The partial decay widths of the heavy Majorana neutrinos Ni are given by 

r tW L = r(Ar . ^ g - w +) = T{Ni i+W^ = _^_\ Ve .\2 M 3 {1 _ /U . Ty) 2 ) (24) 

T iW T _ r(jv . ^ rw -+) = ^|y £ .|2 M . (1 _ Ai . w)2) (25) 

OZ7T 

r „,Z £ = r( ^ ^ nZL) = ^!^|y £ .|2 M 3 (1 _ )2) (26) 

64vrA% 

r , f z T = r(jV . ^ ^ = ^!^|y £ .|2 M . (1 _ )2 (2?) 

32vrc^ 

where /iy = Mj/Mf. If is heavier than the Higgs bosons /i and (see Appendix C for the properties 
of the Higgs bosons in the B — L extension of the SM), one has the additional channels 

V -ih _ T{N . ^ VgK) = 9 2 \v^ M f{\ - fi ih ) 2 cos 2 O , (28) 

64vrM^ 

r „,ff = r(iVi ^ ^ = 9 \ Va fM?(l - vw) 2 sin 2 d . (29) 

64vrM^ 

Therefore, the total width for A^ is given by 

Y N _ ^2 ^2T £Wl + 2T £Wt + r^ Zi + T^ Zt + T^ + T UeH ^j (30) 
At a high mass of Mjy, the branching ratios of the leading channels go like 

T(rw£) ~ T(£ + W£) » r(i/Zt) » r(z//t + i/JT). (31) 

As discussed above, the lepton-flavor contents of N decays will be different in each neutrino spectrum. 
Here, we also study this issue in great detail for cases I and II. In order to search for the events with best 
reconstruction, we will only consider the A" decay to charged leptons plus a W^. 



1. Decays in Case I: Degenerate Heavy Neutrinos 

In Fig. [5] we show the impact of the neutrino masses and mixing angles on the branching fractions of the 
sum of the degenerate neutrinos N% (i = 1, 2, 3) decaying into e, /i, r lepton plus W boson, respectively, 
with the left panels for the Normal Hierarchy (NH) and the right panels of the Inverted Hierarchy (IH), 
assuming vanishing Majorana phases. Qualitatively, it follows the relations in Eq. (fT4l 

BR{^ ± W t ),BR{t ± W t )^ BR(e ± W T ) for NH, 

BRie^W^) > BR^W^^BR^W^) for IH. (32) 
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M N =300 GeV, M„=120 GeV 
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FIG. 5: Branching fractions of degenerate neutrinos iVj — ► £ + W + £ W + (£ = e, /U, r) for NH and IH versus 
lightest neutrino mass with Mjv = 300 GeV and = 120 GeV, assuming vanishing Majorana phases. 



The branching fraction can differ by one order of magnitude in NH case; and about a factor of few in the 
IH spectrum. As one expects that all these channels are quite similar when the neutrino spectrum is quasi- 
degenerate, mi ~ mi ~ nT-3 > 0.05 eV. Therefore, in this simple case one can hope that if the heavy 
neutrino decays are observed in future experiments one should be able to distinguish the neutrino spectrum. 



2. Decays in Case II: Non-degenerate Heavy Neutrinos 

For non-degenerate neutrino spectra we once again study the simple choice: Case Ha U = I. We show 
the branching fractions of processes TVj — > £ + W~ + £~W + {£ = e,/i, r, i = 1,2,3) corresponding to 
the lightest neutrino mass for NH and IH for Mj = 300 GeV in Fig. [6] As noted earlier, in this simplest 
case all \Va\ 2 (£ = e,/j,,r) are proportional to mj. Therefore the branching ratio of iVj — ► £^W^ for 
each lepton flavor is independent of neutrino mass and thus universal for both NH and IH. Although we 
cannot distinguish the neutrino mass hierarchy, we still can tell the difference of the three heavy Majorana 
neutrinos according to different SM lepton flavors in final states of their dominant decay channels. One has 

Si?(e ± W T ) > BR^W^), BR^W^) for N 1: 
BR{e ± W T ) « BR{^W T ) « BR^W*) for N 2 , 
BRfji^W*), BR^W*) > 5i?(e ± VF T ) for N 3 . 
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FIG. 6: Branching fractions of process Ni — ► £ + W + i W + (£ = e, (i,r,i = 1, 2, 3) versus the lightest neutrino 
mass for NH and IH in Case Ha (fl = I), when M t = 300 GeV and M h = 120 GcV. 

This follows closely to the mixing strengths of the light neutrinos in the previous section. 

As discussed previously, Case lib Q, = I g is identical to the above if we identify Ni «-> N%. More 
involved case for O may be some form of superposition of the three decay patterns, that is to be tested 
experimentally by the flavor combinations. 
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FIG. 7: Branching fractions of N± — > £ + W +£ W + versus Majorana phase $2 for NH and $1 for IH in the general 
non-degenerate case when M\ = 300 GcV, M/j = 120 GeV, with random selection of the £1 matrix elements. 

B. Impact of Majorana Phases in Heavy Majorana Neutrino Decays 

In our previous discussion we have shown that the mixings | Vin I 2 are independent of Majorana phases 
in both Case Ha and lib (as well for an Q, with unity as entries). In general, the iV, decay rates depend on 
only one Majorana phase $2 ($1) when mu^\ « and S13 = in the NH (IH) case as shown explicitly 
in the appendices. In Fig. [7J we show the dependence of N\ decay branching fractions for general non- 
degenerate case on Majorana phases $2 and <J>i in NH and IH, with random selection of the Q matrix 
elements. The dependence of N2 and N3 decays on Majorana phases are almost the same as that of Ni. 
The branching fractions of fi^W^, t^W^ (e^W^) are typically dominant over all the range of $2 ($1) 
in NH (IH). The dependence on the phases for the leading channels are rather weak and it is thus hard to 
extract the phase information from heavy Majorana neutrino decay. Some typical situations may be similar 
to the cases discussed in Ref. J17I1 . and we will not pursue further for the phase effects. 



C. Total Decay Width of Heavy Majorana Neutrino 

To complete this section about the heavy Majorana neutrino properties, we study their total decay widths, 
which are proportional to M„M^/M^. In Fig.[U we plot the total width (left axis) and decay length (right 
axis) for N versus under the general non-degenerate case with random selection of the O matrix 
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FIG. 8: The total width and decay length of N in the general non-degenerate case, when the lightest neutrino mass 
10~ 4 eV < mi(3) < 0.4 eV, Mh = 120 GcV and fl = R12R13R23 with random selection of the matrix elements 
and . 

elements (similar for NH and IH). There is a large spread for the possible ranges of the decay lengths, 
governed by the mixing parameters. Although not generally considered as long-lived for large mass, the N 
decay lengths may be typically in the range of /jm-cm, and their decays could lead to a visible displaced 
vertex in the detector at the LHC. 

When considering a specific model-parameter setting, we plot the total width (left axis) and decay length 
(right axis) in Fig.|9j for JVj versus M N for M h = 120 GeV in NH and IH under Case Ha with tt = I. One 
of the generic features for all iVj and both NH and IH is a typical lower limit for their lifetime (or decay 
length). For instance, the typical decay length for Mjv > 600 GeV is above 1 fim. For smaller values of 
M/v, the heavy Majorana neutrinos can be long-lived in the detector scale, making the signatures detectable 
at the secondary vertex. In fact, this feature remains in a majority part of the parameter space. In particular, 
because in this case all \Vei\ 2 are proportional to mi, the lifetimes of N\ in NH and N3 in IH could be 
infinite when neglecting the lightest neutrino mass in whole Majorana neutrino mass range. It is interesting 
to note that there is a clear difference between the NH and IH scenarios: the lifetime of Ni in IH and JV 3 
in NH has a narrowly predicted range within one order of magnitude, about 10 /im for M/y = 400 GeV. If 
this is indeed observed, it could serve as an indication to distinguish the models. The lifetimes of N2 in NH 
and IH are almost the same. 

For Case lib with an off-diagonal matrix, the lifetime features of N\ and N% are also interchanged 
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FIG. 9: Total width and decay length of Ni(i = 1, 2, 3) for NH and IH in Case Ha (ft = I), when the lightest neutrino 

mass 1CT 4 eV < m 1(3) < 0.4 eV and M h = 120 GeV. 



with each other and those of N2 are still the same as Case Ha. 
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IV. HEAVY MAJORANA NEUTRINOS AND THE TEST OF TYPE I SEESAW AT THE LHC 

In order to study the prediction for the lepton flavor correlations with heavy Majorana neutrino and its 
LNV decay processes, the ideal production channels are the Drell-Yan processes via SM gauge bosons, 
pp — > W — > N£, pp — > Z — > NN. However, the gauge couplings to N are highly suppressed to the order 
0{m v /M^) |8J]. The situation is very different in the case of the minimal B — L extension of the SM (see 
appendix C) where one can produce the heavy neutrinos through the Z' in the theory. 

A. Gauge Boson Properties: Z' 

In the limit where there is no mixing between the two Abelian sector of the minimal B — L extension of 
the SM (see appendix C, e = in Eq. (IC3I )). the mass of the new gauge boson Z' is given by 

M z > = 2g BL v s . (33) 

To satisfy the experimental lower bound, Mz'/gBL > (5 — 10) TeV, it is sufficient to assume that v$ > 
2.5 — 5 TeV. The relevant interactions to matter are given by 

9bl % (Q% L [u-fu + drfd] + Qbl [e^e + v L ^v L + vr^vr]) , (34) 

where the B — L charges are assigned to be Q q BL = 1/3, and Q BL = — 1. 



There has been a lot of work on the heavy neutral gauge bosons. For a recent review, see Ref. 11811 . 
and recent studies of Z' at the Tevatron and LHC 11%. For a recent consideration of the phenomenological 
aspects of the B — L model, see |9|1. The expressions for the possible decays of the Z' are given by 

Y{Z> -> //) = 9 BL ^C f (Q f BL ) 2 (l + 2^pj fj f , (35) 

r(Z'^J2 v m"m) = Zg 2 BL ^C v (Q BL )\ (36) 

m 

T(Z' -» N m N m ) = gl L ^-C N (Q BL f (3%. (37) 



where f = £,q, the couplings Q ;I/j jv = 1, C q = 3, and = y 1 — Amj/M^, is the speed of particle i. 
Note that the decay width to Majorana particles is of a threshold behavior /3 3 , and is half of that for a Dirac 
particle. Well above the threshold, the Z' decay branching fractions take the simple ratios for the final states 

e,fi,r u...t li2,3 

y™;/™ : N-i N-i =3:2: 

2 2 
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FIG. 10: Total decay width of Z' , when v s = 3 TeV. 



We show in Fig.[[0]the results for the case vs = 3 TeV. It scales as 



T(tot) « 0.2 g% L M z > < 0.05 



v s 



Mz>. 



(39) 



Notice that this Z' has the property that its coupling to quarks is suppressed with respect to the couplings to 
leptons. As it well known the Z' in Left-Right symmetric theories has different properties from the B — L 
case studied here. Then, from the standard analysis where one uses the leptonic channels and the channels 
into heavy quarks one can distinguish the B — L case from the rest easily. 

It is important to emphasize that in the case of the B — L SM, one gets an upper bound on the mass of 
the heavy neutrinos Mn < Mz> / (2\/29_bl) (see appendix C for details). 

B. Heavy Majorana Neutrino Production through Z' mediation at the LHC 

We are interested in the production of two heavy neutrinos. Since in this model one has a dynamical 
mechanism for B — L breaking, there is a production mechanism through the Z'. Then, we are interested 
in the mechanism 



pp^ Z' -> N\N\. 



(40) 
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FIG. 11: Heavy Majorana neutrino pair production total cross section at the LHC versus its mass. The solid, dashed 
and dotted curves are for Mz> = 1,1.5,2 TeV respectively, when vs = 3 TeV, (a) for U (1)b-l coupling and (b) for 
U(l)x coupling, as given in Tables HI and HI1 

The partem level cross section for this process is 

dt ~^N C 9 ( S -M 2 ,) 2 + M|,r|, [{t Mn} +{U Mn) 2sMn1 

(41) 

where t = (p q — pn) 2 - The total cross section versus heavy Majorana neutrino mass at the LHC is plotted 
in Fig.QTJ assuming v$ = 3 TeV with (a) for U(1)b-l coupling and (b) for U(l)x coupling, as given in 
Tables HI and III1 We see that the production cross sections are quite sizable, typically of the order of 10 — 100 
fb. The cross section drops sharply after reaching the kinematical threshold 2Mjy > Mz> ■ 
The Majorana signals for AL = 2 decay of N\ are 

NxNx -> l ± £ ± W T W T , £ = e,fi,T (42) 

To confirm the important feature of lepton number violation, we demand the W's decay hadronically. The 
overall branching fraction to be included becomes 

BR(JViM ^4 jets) « 2 • (i) 2 • (^) 2 = -L. (43) 

Note that there are also accompanying clean channels like l~^t^ + 4 jets, that are not lepton-number 
violating and we do not include for the rest of the analysis. 
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We would like to reiterate that in a significant range of the parameter space of Mn and mixings, the N 
decay could lead to distinctive signatures with a decay length longer than 10 fim, resulting in secondary dis- 
placed vertices. This may yield essentially background-free signal for iV's. Nevertheless, we now explore 
the signal observability according to the different lepton flavors without relying on the displaced vertex 
considerations. 

For our numerical analyses, we adopt the CTEQ6L1 partem distribution function if^oll . We evaluate the 



SM backgrounds by using the automatic package Madgraph [21Q . We work in the parton-level, but simulate 



the detector effects by the kinematical acceptance and employ the Gaussian smearing for the electromag- 



netic and hadronic energies 12211 



C. iViiVi ^■£ ± ^ ± + 4jets(^ = e,/i) 
We start from the cleanest channels with e, u in the final state from N\ decay. We employ the following 



basic acceptance cuts for the event selection j22ll 



pr(£) > 15 GeV, \rj(£)\ < 2.5, (44) 
Pr(j) > 25 GeV, \t](j)\ < 3.0, (45) 
AR j:j > 0.3, AR je , AR U > 0.4 (46) 

The rather loose cuts on the separations AR are designed to keep the signal events for a heavier Z' and a 
lighter N which is fast moving and thus yields collimated decay products of a lepton and two jets. We plot 
the minimal isolation ARJj in of two jets and AR™ n of one jet and one charged lepton for Mz> = 1 TeV 
and Mjy = 100, 200 GeV, respectively, in Fig.[l2] One can see that for M N > 200 GeV with M z > = 1 TeV 
the signal consists of well-isolated one pair of same-sign leptons of arbitrary e, /i flavor combinations plus 
four light jets. 

To simulate the detector effects on the energy-momentum measurements, we smear the electromagnetic 
energy, the electromagnetic energy and jet energy by a Gaussian distribution whose width is parameterized 



as |22[] 



^ = J^g— b cal , a cal = 10%, b cal = 0.7%, (47) 
& y/E/GeV 

= ^ b had , a had = 50%, b had = 3%. (48) 

E y'E/GeV 

In principle, there is no genuine SM background to the lepton-number violating processes. The leading 
SM background to our signal is from decays of two like-sign W's to leptons. For instance, the leading 



21 




0.5 1 1.5 2 0.5 1 1.5 2 

AR min AR min 

FIG. 12: AR™ m and Ai?™ in for M Nl = 100 GcV (left) and M Nl = 200 GcV (right), with M z > = 1 TcV. 

reducible background to our signal is 

pp -> tiW ± -» PF ± W ± ji&6. (49) 

The QCD processes jjjjW ± W ± ,jjW ± W ± W Zfl are much smaller. This is estimated based on the fact 
that QCD jjW^W^ — > jji^i^J^x is about 15 fb. With an additional and 6 body phase space or one 
more suppression, they are much smaller than tiW^. Other EW backgrounds WWWW, WWW Z are 
also neglectable. Although the background rates are large to begin with, the kinematics is quite different 
between the signal and the backgrounds. We outline the characteristics and propose some judicious cuts as 
follows. 

• The SM backgrounds always come with W pair decays with missing neutrinos. To suppress back- 
grounds, we veto the events with large missing energy < 20 GeV. 

• We choose the two pairs with nearly equal masses from the six dijet combinations as the two hadronic 
W's and take W boson reconstruction as \Mjj — Myy \ < 15 GeV. The efficiency is very high. 

• In order to select the correct lepton and two jets combination and reconstruct Ni , we take advantage 
of the feature that the two heavy neutrinos have equal masses Mi 1 j 1 j 2 = Mc 2 j. s j i . In practice, we 
impose \Me 1 j 1 j 2 — Mg 2 j 3 j A \ < M^/25. This helps for the background reduction. 

The production cross section of N%N% signal with the basic cuts (solid curve) and all of the cuts above 
(dashed curve) are plotted in Fig. \Y3\ where branching fractions for Ni decay to charged leptons are not 
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FIG. 13: Production cross section of N\N\ with basic cuts and hard final states cuts. Branching fractions for N± 
decay to charged leptons are not included; while W decays to 2 jets are included. For comparison, the background 
process is also included with the sequential cuts as indicated. 

included; while W decays to 2 jets are included. For comparison, the background process of tiW^ is also 
included with the sequential cuts as indicated. The background is suppressed substantially. 

When performing the signal significance analysis, we look for the resonance in the mass distributions 
of £jj and 2£4j. If we look at mass window of \M£ 1 j 1 j 2 ^ 2 j 3 j 4 — Mjv x | < and \M 2 e+4j — Mz> \ < 

Mz' /SO, the background will be at a negligible level. 

D. iViJVi -*-r ± £ ± +4jets 

The previous section sets the stage for the analyses in the following sections. Most of the issues for event 
selection and detector acceptance will remain the same for the following studies. The next presentations 
will thus be sketchy and mainly outlining the new features, in particular the r-reconstruction and the mass 
resonances. 

The r lepton final state from heavy Majorana neutrino decay plays an important role in distinguishing 
different neutrino mass patterns. Its identification and reconstruction are different from e, \x final states 
because a r decays promptly and there will always be missing neutrinos in r decay products. In practice 
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when selecting events with r's, we require a minimal missing transverse energy 

Xt > 20 GeV. (50) 

This will effectively separate them from the lljjjj type of signal events. 

We first note that all the r's are very energetic from the decay of a few hundred GeV Ni. The missing 
momentum will be along the direction of the charged track. We thus assume the momentum of the missing 
neutrinos to be reconstructed by 

"^(invisible) = Kp (track). (51) 
Identifying pr (invisible) with the measured JZt, we thus obtain the r momentum by 

~Pt(t) = + Xt, P L (t) = p L (i) + ^PlW- 

The Ni pair kinematics is thus fully reconstructed. The reconstructed invariant masses of M(£jj) and 
M(rjj) are plotted in Fig. |T4j We see that M(rjj) distribution (dotted curve) is slightly broader as antici- 
pated. The rather narrow mass peak of the Ijj system nevertheless serves as the most distinctive kinematical 
feature for the signal identification. Invariant masses of M{rl + 4j) are also plotted in Fig.O Although 
the existence of missing energy in the signal makes the background separation more involved, the resonant 
mass reconstruction proves to be high efficient and the backgrounds can still be suppressed to a negligible 
level. 

E. -> r ± r ± + 4 jets 

For rrjjjj events with two r's, we generalize the momenta reconstruction to 

"^(invisible) = n± p (tracki) + k 2 ~p > (track 2 ). (52) 

The proportionality constants k\ , k 2 can be determined from the missing energy measurement as long as the 
two charge tracks are linearly independent. The N\ pair kinematics can be once again fully reconstructed. 
The reconstructed invariant masses of M(rjj) is plotted in Fig. [15] The nice mass peaks of the rjj system 
at Mn 1 and the rrjj system at Mz> make the signal stand out of the SM backgrounds. 

It is important to note a difference between the leptons from the primary N\ decay and from the r decay: 
the latter is much softer. In Fig. [16] we show the px distribution of the softer lepton from the N\ and r decays 
in the events of lljjjj, li~jjjj and rrjjjj. This feature could provide additional discrimination power to 
separate the three different leptonic channels if needed to fit the flavor structure for a underlying theory. 
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FIG. 14: Reconstructed invariant mass of M(jj£), M(jjr) for Mjv x = 300, 600 GeV, respectively (left), and 
M(rUj) fox M z > = 1, 1.5 TcV (right). 

F. Measuring Branching Fractions and Probing the Neutrino Mass Patterns 

So far, we have only studied the characteristic features of the signal and backgrounds for the leading 
channels and have not included the proper branching fractions for the individual lepton flavors. For illustra- 
tion, consider first the cleanest channel, N\N\ — ► e e^jjjj. The number of events is written as 

N = L x a(pp -> NxNx) x 2 BR 2 (iVi -» e + J^-)(-) 2 , (53) 

9 

where L is the integrated luminosity and the factor (6/9) is due the the W hadronic decay. Given a sufficient 
number of events N, the mass of N\ is determined by the invariant mass of lepton and jet Mgjj. We thus 
predict the corresponding production rate a(pp — ► NiN\) for this given mass. The only unknown in the 
Eq. (|53l ) is the decay branching fraction. 

We present the event contours in the BR— M^r plane in Fig. [F7]fbr 100 fb _1 luminosity and degenerate 
case with (a) Mz> = 1 TeV and (b) Mz> = 1.5 TeV including all the judicious cuts described earlier, with 
which the backgrounds are insignificant. 

In Fig. [l7](c) and (d), we show the event contours in the BR— plane, for 100 fb" 1 luminosity and 
non-degenerate case including all the judicious cuts described earlier. We see that the reach to a low BR can 
be quite encouraging. 

As we presented earlier, the N\ decay branching fractions and the light neutrino mass matrix are directly 
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FIG. 15: Reconstructed invariant mass of M(jjr) for M^ 1 = 300, 600 GeV, respectively (left), and M(rr4j) for 
M z > = 1, 1.5 TcV (right). 



correlated. Measuring the BR's of different flavor combinations becomes crucial in understanding the neu- 
trino mass pattern and thus the mass generation mechanism. In the degenerate case, we have the prediction 
for the flavor combinations 



BR(NN -> UWW) 



2 x (23%) 2 for NH : (/x ± + r ± )(^ ± + t±)WW, 
2 x (13%) 2 for IH : e ± e ± WW, (54) 
2 x (17%) 2 for QD : (e ± + n± + r ± )(e ± + /i± + r^WW, 

for $i = $2 = 0, independent of the matrix fl. On the other hand, for the non-degenerate situation, the 
flavor prediction is like 

' 2 x (20%) 2 for Ni : e ± e ± WW ) 
2x(17%) 2 fovN 2 :(e ± +fi ± + T ± ){e ± +fi ± + T ± )WW, (55) 
2 x (23%) 2 for N 3 : (/x ± + r ± )(^ ± + r^WW, 

for fi = /, independent of the neutrino mass patterns as well as <&i, $2. These predictions are the conse- 
quence from the low energy oscillation experiments and this are subject to test at the LHC to confirm the 
theory. 

In Fig. [18] we show the event contours in the Mz> — Mjy plane, for (a) production of ^ iVj in NH (solid 
curve), IH (dashed curve) and QD (dotted curve) for the degenerate case and (b) production of N\ (solid 



BR(iV7V UWW) 
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FIG. 16: pt distribution of the softer lepton from the N\ and r decays in the events of Mjjjj, Irjjjj and rrjjjj, 
for a N x mass 300 GeV and Z' mass 1.5 TeV. 

curve), N2 (dashed curve) and N3 (dotted curve) for the non-degenerate case with 100 fb _1 luminosity, 10 
events numbers and branching fractions of heavy neutrinos predicted in Eqs. (l54l) and (I55T ). The values of 
Mz> and on the left-hand side of the curves would give more than 10 events for 100 fb" 1 luminosity 
and more accessible heavy neutrino decay branching fractions at the LHC. 
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FIG. 17: Event contours in the BR— Mn plane at the LHC with an integrated luminosity 100 fb 1 for degenerate 
case £\=i 2 3 NiN t -> e + e + W~W~ with (a) M z > = 1 TeV, (b) M z > = 1.5 TeV, and for non-degenerate case 
A^iA^i — > e + e + W~W~ with (c) Mz< = 1 TeV and (d) M^/ = 1.5 TeV, including all the judicious cuts describe in 
the early sections. 

V. SUMMARY 

In this article we have investigated the possibility to test the so-called Type I seesaw mechanism for 
neutrino masses at the CERN Large Hadron Collider in the context of two simple extensions of the Standard 
Model where B — L is part of the gauge symmetry. We have studied in great detail the predictions of the 
right-handed neutrino decays in each spectrum for neutrino masses showing the most optimistic scenarios 
where one could hope to distinguish the spectrum using the properties of the decays. 

We have found the following interesting results: 
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FIG. 18: Event contours in the Mz 1 — Mjv plane at the LHC with an integrated luminosity 100 fb _1 and 10 events 
number for (a) degenerate case $^ i=1 2 3 NiNi — > l^l^WW for NH, IH, and QD, and (b) non-degenerate case for 
Ni,N 2 and N 3 , all as predicted in Eqs. (0) and dS). 

• Working in the context of two simple extensions of the Standard Model with a local gauge symmetry 
B — L or X = Y — j(B — L), one can produce the heavy neutrinos through the Z' gauge boson 
in each scenario. In both cases one has a dynamical mechanism for the generation of heavy neutrino 
masses, related to Mz>. 

• In the case where the heavy neutrinos are degenerate, we show the possibility to distinguish the 
neutrino spectrum. The branching fractions can differ by one order of magnitude in NH case 
with BR(h ± W t ),BR(t ± W t ) » BR{e ± W T ), and a factor of a few in the IH spectrum with 
BR{e ± W T ) > BR(^L ± W Zf ), 5i2(r ± W =F ) when the Majorana phases are ignored. As one expects, 
all these channels are quite similar when the neutrino spectrum is quasi-degenerate, mi ~ m-2 ~ 
m 3 > 0.05 eV. 

• In the case when fi is an identity matrix or with only unity entries generally, we find: BRfe^W^) > 
BRfjj^W*), BRfj^W*) for Ni decay, BR^W^) w BRfji^W*) ~ BR^W*) for N 2 and 
BRQjfiW*), BR^WT) > J B J R(e ± W =F ) for N 3 in both NH and IH. The branching fractions in 
these cases are independent of Majorana phases. 

• In general, the form of 17 governs heavy neutrino decay patterns. Future tests on the flavor combina- 
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tions of SM charged leptons would reveal the specific model structure. 

• The above-studied AL = 2 channels can take the search to Mz' ~ 2 TeV at the LHC. The sensitivity 
to the leptonic branching fractions of N decay can be about 10%. 

• In particular, in a significant part of the parameter space of M n and the mixings, the N decay could 
lead to distinctive signatures with secondary displaced vertices, yielding essentially background-free 
signal for N's. 
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APPENDIX A: NEUTRINO MASSES AND MIXINGS 



The Type I seesaw scheme introduces right-handed neutrino states, in addition to the SM matter contents. 
For definitiveness, we add three of them, v l R (i = 1, 2, 3), which are said to be sterile since they do not carry 
any SM gauge quantum numbers. The SM gauge invariant and renormalizable interactions to generate 
neutrino masses include both Dirac as well as Majorana terms 



- Ci = l L Y v u Hu R + -(u c ) L M N v R + h.c. (Al) 

where Y^, M/v are 3 x 3 matrices in the generation space, with H = ia^H* and H T = (H + H°Y Once 

H gets the vacuum expectation value (H) = vo/y/2, the neutrinos acquire Dirac masses m D = Y® vq/V2, 

~£™=2 { Vl mD VR + m D (" C )r + W)l M n vr) + h-c. (A2) 
To diagonalize the mass matrix for neutrinos we introduce a 6 x 6 unitary transformation 

N \ Uh \ , N = I V 

S» c )l ) \ (v c )l \ V C U C 

Then, 

N t I " "' u J N* = I " I , (A4) 






m D l 


I N* = I 


( m \ 






V m D 


M N J 




I M I 
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or explicitly, 

V^m T D U* + U ] m D Vc + V&M N Vc = m, (A5) 

U ] c m T D V* + V^m D Uc + U ] c M N Uc = M, (A6) 

V^m^F* + Wm D Uc + V^Mivt/o = °> ( A7 ) 

where m = diag(mi, rri2, m^j and M = diag(M\, M2, M3) are diagonal matrices of the mass eigenvalues. 
In the limit Mjv > roc we have 

2 

Then we have three light neutrinos and three heavy ones, all Majorana-type. Note that U and V mix the 
light and heavy neutrinos, respectively, into the active weak interaction eigenstates of neutrinos. The mixing 
elements are typically like 

U 2 ^0(l), V 2 ^^, (A9) 

and the unitarity conditions read 

UU ] + VV* = U ] U + VXV C = V C VX + UcUl = V ] V + Ulu c = I, (A10) 



r V c l 



UVX + VUt = U^V + VXU C = 0. (Al 1) 



It can be shown that the following relations hold 

V^m T D - mU T = 0, m D Uc - VM = 0. (A12) 
Assuming that a 3 x 3 matrix E diagonalizes the mass matrix of the charged leptons, we then define 

E ] U = Vpmns, E^V = V iN , (A13) 

VpmNsVpmNS + V tN V lN = 7 ' ( A14 ) 

where Vp m ns and Vepj describe the transitions between the light neutrino and heavy neutrino to the charged 
leptons, respectively, via the weak charged currents. Note that the identification of Vpmns to the PMNS 
matrix is only approximate. We then obtain an important relation among the physical quantities 

V e * N M V} N = -Vhms m V PMNS . (A15) 

Although the masses and mixings of the light neutrinos on the right-handed side can be measured from the 
oscillation experiments, it is quite involved to solve for Vex via this set of quadratic equations. Absorbing 
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the minus sign on the right side of the above equation in the definition of Vgjy one can write down a formal 
solution with the help of an auxiliary matrix 



V eN = VpMNS^VLM- 1 / 2 , 
where Q is an orthogonal complex matrix which can be parameterized as 

0(^21,1031, ^32) = #12(1021) #13(1031) #23(1032), 

with 



i?12 



U21 —W21 ^ 
w 2 i u 2 i 
1 



Ri 



1 u 31 -to 3 i ^ 

1 
\ W31 u 3 i 1 



R 2 



I 1 ^ 

« 32 -U>32 
y W 32 "32 J 



where Uij = — w 2 j and — 1 < Wy < 1 when the matrix f2 is real. 



(A16) 



(A17) 



(A18) 



APPENDIX B: EXPLICIT EXPRESSIONS FOR THE MIXINGS Vw 



1. Case I: Degenerate Heavy Neutrinos 



From Eq. ( IA15I ). assuming degenerate heavy neutrinos, we have 

N2_,{V eN ) = c 13 s 12 m 2 + c 12 c 13 e 1 m 1 + s 13 e ( 2> m 3 , 

N 

^n^{V* n ) 2 = (C12C23 - si2Si3S23e~* <5 ) 2 m 2 + (c 23 si2 + ci2Si 3 S23e" 8,5 ) 2 e" l * 1 TOi 



JV 



2 2 -i* 2 

c 13 s 23 e 7713, 



M N J2(V; N ) 2 = (ci 2 s 23 + c 23 si2Si 3 e lS ) 2 m 2 + (s 12 s 23 - c 12 C2 3 s 13 e lS fe l$1 



771l 



c 13 c 23 e 4<1>2m 3- 



(Bl) 



(B2) 



(B3) 



2. Case II: Non-degenerate Heavy Neutrinos 



The general expressions for the mixing between the charged leptons and heavy neutrinos, in terms of the neutrino 
oscillation parameters and the unknown matrix Q, are given by 

V e iy/Mi = V^c 13Sl2 W2i^/l - w 2 31 + V^rc 12 c 13V /(l - w 2 21 )(l - w 2 31 )e^' 2 + V^ S i3«;3ie l( * 2/2 - 5) , 

(B4) 



V^iV-Mi = y/m^(ci2C2 3 - si 2 s 13 s 23 e l )w 2 i J 1 - w 31 



+ ^(-^2^3 - c 12 s 13 S2 3 e i5 )^l - w 2 21 ){\ - w 2 3i y^/ 2 + 



/m 3 ci 3 s 23 w 31 e 



i* 2 /2 



V T i^fM~i = ^/m2~(-c 12 s 23 - s 12 s 13 C2 3 e lS )w2iJl - w 31 



+ V^T(si2S 23 ~ c 12 .s 13 c 23 e i;5 ) v /(l - w 2 21 ){l ~ w 2 x y^l 2 + V^c^^i e 4 * 2/2 . 



(B5) 



(B6) 
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V f WM 2 = 



V tt2 x/M 2 



+ 



frn2~{c 12 c 23 - s 12 s 13 s 23 e l5 ){-w 21 w 31 w 3 2 + \J {I - - wf 2 )) 

/ TOi(-si 2 c 2 3 - ci 2 si 3 S23e 4 ' 5 )(-W32W3iy / l - «4l - w 2 i^J 1 - wf 2 )e l * l/2 



/ to 3 "ci 3 S23W 3 2\/ 1 - ui^e 4 * 2 / 2 , 



= \An2(-Cl2S23 - Sl2Sl 3 C 23 e l5 )(~W 2 lW 3 lW 3 2 + \J (1 ~ ~ W 32 )) 

+ V™l( S 12 s 23 - Cl2Sl 3 C23e l ' 5 )(-W32'^3iy / l ~ W%1 - W 2 lJl - W^e^ 1 ^ 



+ y/rmcuc^w^J 1 - w 31 e' l ' s>2/2 



(B7) 



(B8) 



(B9) 



VWMg 



/ TO2Ci3Si 2 (-W32y / l - - W2lW3lJl ~ W32) 

/ mici 2 ci 3 {w 21 w 32 - w 3 i\J (1 - w^^l - wl 2 ))e 4 * l/2 

^3^(1 - <)(i - t^&K^/ 2 - 5 ), 

/ m2(ci2C 2 3 - si2Si3S 23 e l ' 5 )(-w 3 2y / l - w 21 - lV 2 iW 3 i^Jl - wj 2 ) 
r rn\{-s 12 c 23 - c 12 si 3 s 23 e l5 ){w 32 w 2 i - w 31 \j (\ - w\-y){\ - wf 2 ))e l * l/2 
^c 13 s 23 ^/(l - w 2 31 )(l - ^ 2 2 )e 4 * 2/2 , (Bll) 



(BIO) 



V t3 \JM 3 = y/m2~(-ci2S 23 - s 12 s 13 C2 3 e l ){~w 32 \J 1 - w\ x - w 2 iw 31 y 1 - w 2 2 ) 

+ \Ml(si2S 2 3 - Ci 2 S 13 C 23 e l5 )(w 32 W 2 l ~ W 3 i\J(l ~ w\^)(\ ~ «7 2 2 ))e I * l/2 

+ 0^c 13 C2 3V /(l - w 2 31 )(l - w 2 2 )e^-/ 2 . (B12) 
We now present the two cases according to the light neutrino mass spectra, assuming mjp) sa and s 13 = 0. 
• Normal Hierarchy: 

Under the good approximations mi ps and s i3 = 0, one finds the following expressions 



M x \V el \ 2 « JAm^^xa^iVl-^) 2 , (B13) 



Mi|V M i| 2 « IV Am 2iCi2C23^2iVl-i«ii + (/ Am 2 ! s 23 w 3 i e 1 * 2 /2 1 2 , (B14) 



^ilKil 2 « |^Am§ 1 d 2 *23«Wl- tu 3i- \7Am 2 lS2 3W3ie 8 * 2/2 | 2 . (B15) 



M 2 \V e2 \ 2 



M 2 \V T2 \ 2 



\J Am^sl 2 (-W2lW S lW 3 2 + - w 2 21 ){l - w 2 2 )f 



M 2 |y M2 | 2 « |yA^Cl 2 C2 3 (-K;2lW31^32 + ^(1 - - W§ 2 )) 



| $J Amj 1 c 12 s 23 (~w 21 w 31 w 32 + \J '(1 - i«2i)(l - ^32)) 



M 3 |y e3 12 



I 2 ~ \J ^m 2 2l s 2 l2 {~w 32 ^Jl - w\\ - w 2 iw 31 ^l- w 2 2 ) 2 



M 3 |V m3 | 5 



M 3 \V T3 \ 2 



Am? 21 c 12 c 23 (~w 32 ^J 1 - - w 2 iw 3 i^Jl - u>l 2 ) 
tf&mgs^l - wl)(l - w 2 2 )e^ 2 \ 2 



Am 2 21 c 12 s 23 (~w 32 \/ 1 - u- 2 ^ - w 2 iw 3 i^Jl - wl 2 ) 
yA^ C23V /(l - ^(1 - W3 2 2 )e^ /2 | 2 . 



Inverted Hierarchy: 

Under the approximations m 3 ps and si 3 = 0, we have 



|Am3 1 |si 2 W2i\/l 



+ ^Aml^il-wl^l-w^y^/'l 2 , 



MxlKil 2 



M^V^ 1 « |^/A77i2 1 + |Aml 1 |ci 2 c 23 W2iyi-wii 

- ^A^4| S1 2C2 3) /(1 - a&Xl - u&K*^ 2 ! 2 , 
MxlKil 2 « 



M2IK2I 2 



m^j + | Am 31 |ci 2 s 23 W2i yl-w^ 
^\A m 2 31 \ Sl2 s 23 ^{l - w 2 21 )(l - w^e^ 2 ] 2 . 

\J Am^ + |Am| 1 |si 2 (-w 2 iw 3 iw 3 2 + \J (I - - w| 2 )) 

+ ^A^c 12 (- W32 ^/l - < - w 2 x/l - w 2 2 )e^/ 2 \ 2 , 
M 2 |T^ 2 | 2 « l^/Am 2 ! + |A77T.§ 1 |c 12 c 23 (-W2iW3iW32 + ^(1 - tu£0(l - w 2 2 )) 
- \/|Am^ 1 |si 2 c 23 (-w 3 iu; 3 2y / l - loli - w 2 iy / l - w| 2 )e ! * l/2 | 2 , 
I 2 ~ I ^/ Am^ + |Ato| 1 |ci 2 s 2 3(-W2iW3iw 3 2 + ^(1 - - wf 2 )) 

|Am| 1 |si 2 s 23 (~u; 3 iw 32 y / l - w\ x - w 21 ^\ - w| 2 )e l * l/2 | 2 . 



M 2 |y r2 |2 
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M 3 \V e3 \ 2 w l^/Am^ + \Aml 1 \s 12 (-w 3 2y/l ~w 2 i -w 21 w 3 i^l - W32) 

+ ^\A^\c 12 (w 21 w 32 - w 31s J(l - w 2 21 )(l - w 2 2 ))e^/ 2 \ 2 , (B28) 
A^IV^I 2 ps \{J Am^ + |Am| 1 |ci2C 23 (-W32y / l - - w 21 w 31 ^l - w 2 2 ) 

V^A^J Sl2 c 23 (^21^32 - ^31 x/C 1 - - ^| 2 ))e 1<I>l/2 | 2 , (B29) 

M3IK3I 2 « I^/Ato^ + \ Aml 1 \c 12 s 23 (-w 32 ^ 1 - w 2 21 - w 21 w 3l ^l- w\ 2 ) 

tf\A^\s 12 s 23 (w 21 w 32 - w 31 ^{l~w 2 21 )(l-w 2 2 )y^ /2 \ 2 . (B30) 



APPENDIX C: U(1) B -l AND C/(l)x EXTENSIONS OF THE STANDARD MODEL 

It is well-known that i? — L is an accidental global symmetry in the Standard Model and its origin is unknown. 
In order to understand the origin of Majorana neutrino masses it is crucial to look for new scenarios where B — L 
can be spontaneously broken. Here we focus on a simple extension of the Standard Model where U(1)b-l is a local 
symmetry and in order to cancel the anomalies one has to introduce three right-handed neutrinos. Therefore, this 



model is based on the gauge symmetry SU{3) C SU{2) L U(1) Y U{\) B -l 
following properties: 



The matter fields have the 




(3, 2, 1/6, 1/3), u R ~ (3, 1, 2/3, 1/3), d R ~ (3,1, -1/3, 1/3), (CI) 



~ (1, 2, -1/2, -1), e R ~ (1, 1, -1, -1), and v R ~ (1, 1,0, -1), (C2) 

L 

where v R are the right-handed neutrinos. Here we use the normalization where Q = T 3 + Y . In order to generate the 
right-handed neutrino masses and break the local B — L symmetry one has to add a new scalar field S ~ (1,1,0, 2). 



1. Interactions and Symmetry Breaking 

In this context the Kinetic terms for the Abelian sector are given by 

'-'gauge — ,r r ptv ^ r r ^ v ^ /«" \^ J J 

where 

= d»B v - d v B», and F 1 ^ = dy0 v - d v B'^. (C4) 

Here B v and B' v are the gauge fields for U(\)y and U(1)b-l, respectively. Since the mixing between the Abelian 
gauge bosons have to be very small we work in the case where e = 0. The Kinetic terms for the matter fields read as: 

C-Kinetic = iQl^D^Ql + iur^D^ur + Mr^ D ^dR 

+ Ul^DJl + ien^D^R + iv R ^D^ R , (C5) 
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where 

D^r = d )l v R - ig B hB'^ R . (C6) 

As we have explained before the Higgs sector is composed of the SM Higgs, H T = (H + . H°), and an extra Higgs, 
S = S R + iSi, which is needed to break B — L. The relevant Lagrangian for the scalar fields is given by 

C-Hi aas = (D^ (D»H) + (D,J3)1 (D»S) - V(H, S) (C7) 

where 

D^S = d^S + i2g BL B^S. (C8) 
The gauge invariant Yukawa interactions of neutrinos are 

yM 

-CX, = Y„ D l L Hv R + ~yVrC vr S + h.c. (C9) 

Once S gets the vacuum expectation value (S) — > vs/V2, B — L is broken and one gets the mass of neutral gauge 
boson Z' = Zb-l with Mz 1 = 2g B LVs from the second kinetic term in Eq. (|C7t , and the mass matrix of right- 
handed neutrino with Mn = Yj^vs / \/2 from Eq. (|C9t . In order not to upset applicability of perturbative theory, we 
require Y^ L < 1 and get an upper bound of the mass of the heavy neutrinos Mjv < Mz> / {2\]~2g B i). 
The scalar potential is given by 

V(H,S) = -m% H^H + X H (H^H) 2 - m| S^S + X s (5 f 5) 2 + a s (H^H) (5 f 5) (CIO) 

where all parameters are real. Notice that this scalar potential has the global symmetry 0(A) B <S> 0(2)s- The mini- 
mization conditions in this case read as 

= v {-m 2 H + Xhv 2 + y v s) > < C11 ) 







v S (-m 2 s + X s v 2 s + ^-vf) . (C12) 



Notice that one can have several vacua but only the case vq ^ and v$ ^ is allowed by the experiment. Now, in 
order to satisfy the condition of minimum one has to satisfy the following condition: 



>^Hasv G + 4X H X s v Q v s + X s a s v s > 0. (C13) 

The potential is bounded from below when XrXs — <x|/4 > 0. Using the minimization conditions above one can find 
the solution in the phenomenological allowed case: 

2 _ 2 (a s m 2 s - 2X s m 2 H ) 



a 2 s - 4X H X s 



2 _ 2 (a s m 2 H - 2X H ™ 



ni 



s) 



> 0, (C14) 



> 0. (C15) 



a 2 - 4X H X S 

Using these conditions one can discuss different cases for the parameters in the Lagrangian. Expressing the numerators 
and the denominator as n\ = as m 2 H — 2 Xh m 2 s , ri2 = as m s ~ 2 Xs m 2 H and d = a| — 4 Xs X B - 

• Imposing v% > one has the case n\ > and d > 0, or n\ < and d < 0. 

• Imposing v 2 > one has n-i > and d > 0, or ri2 < and d < 0. 
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Fields 


Vertices 


Couplings 


Approximations 


Z' 


qmZ' 
Ql = u, q 2 = d 


-iQ q BL9BLj^ 

O q - i 


- 




UZ' 


-iQBL9BLl^ 


- 




I = e, fi, t 


Qbl = -1 






N mi N m2 Z' 


-i{U%U* c - V T V*) mim2 Q" BL g BL r^ 


Hnn^gBL^^ 






-i(WU - V^ c Vc) mi m 2 Q l BL 9BLl^ 


H mi m 2 9BL^^ 








- 




NII-W+ 




- 






2cw mim 2 ' ^ 






T 

^mi-^mj Z 


2c w u m 1 m' 7 ' 1 ^ 






T7lN m h 


• xr D (M» . AC \ 

-iVemPR {-^r c 9o + -^r s e«) 


M N 




ViN m H 


•xt p (Ml ME \ 


-iV lm P R -^s eQ 



TABLE I: Feynman rules for Z' and heavy Majorana neutrino N in SM with U{1)b-l extension, where U = U^V. 

2. Higgs Bosons Properties 

As we have discussed before the Higgs sector of this model is composed of the SM Higgs, H T = (H + , (vq + 
H° + iS,°)/^/2), and an extra Higgs, S = (vs + S° + iSi)/V2, which is needed to break B — L and generate 
neutrino masses. In this context one will have only two CP-even physical Higgses h and H, and the mass matrix for 
the these fields is given by 

M 2 Q = S . (C16) 

y asv v s A s u|/2 - J 

The physical Higgses are defined by 

h \ I cosfti sinfti \ / H° \ 

(C17) 





where the mixing angle is 



COS Vq sm "0 

— sin Oq cos 6*o 



4. oa a s v v s 

tan26» = - * — ^- 7-5 o\. (CI 8) 



It is easy to check that Sj is the Goldstone boson eaten by the Z' in the theory. 



3. Feynman Rules 



We now summarize the Feynman rules for the SM with U(1)b-l ar, d U(l)x extensions in Tables HI and ITfl 
respectively. 
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4. Z' Decays in U(l)x Extension 

The charge of U(l)x is defined as X = Y — 5(B — L)/4 and due to the mixing between U(l)y and U(l)x we 
have the mixing matrix of neutral gauge bosons as below 





1 








\ 



sw cwc' —cws' 
s' d 











) 


{ z>» ) 



(C19) 



where s w (cw) = sin 6\\r (cos 6\y), s'(c') = sin(0')(cos(0')), andtan(2#') = 2g[ yj g\ + g\j (gf + 25g% L ^| — g\ — 
g\). g[ is a free gauge coupling to qualify the mixing between the two U(l) gauge symmetries. One can get their 
mass eigenvalues as 



M A 
M z ,z> 



vo 
2 



gl + al 



i fg? 



9l 



9l 



9i 



sin2^V^+.92 2 



1/2 



(C20) 
(C21) 



The expressions for the possible decays of the Z' are given by 

M z 



r(z'^ff) = 

m 

T{Z'^N m N m ) = 

T(Z' ->W + W-) = 

Y{Z' -> Zh) = 



r(z' -> zh) 



12tt 

n M z > 



Vf\l 



C v {X?f 



2m) 



M% 



A - 



l - 



Am) 
M%, 



1 - 



24tt 
Mz ' r (Y u \< 

aM z ,s' 2 M.%, 
3tan^ WM^r 



'N 



Ml 



Ml 



I H 

Ml 
68M| 



48M 



ir 



M%, 



i M z > r 

48tt Mf 



5Mf 
2 z 



M, 2 



Mf, 



Mf + M* [Ml-MlY 



Ml 



Mf, 



i m Z / r 

48tt M| 



uo^ 2 -^ + v s ce K' 



5Mf 
2 z 



Mf 



A7 2 , 
2 z 



Mf 



Mf, 



Mi, 



4m| 
Mf, 



4M^ 
Mf, 



(M 



(Mf-M f) 2 
Mf, 



(C22) 
(C23) 

(C24) 

(C25) 



(C26) 



(C27) 



where / = .£,?, CIi,„ jJV = 1, C, = 3 and 7/ = + X fR )/2,A f = {-X fh + X fR )/2. 
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